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FREQUENCY  UPSHIFTING  OF  LASER  PULSES  BY 
IONIZATION  FRONTS 


The  propagation  of  electromagnetic  radiation  through  plasmas  is  a  problem  of  general 
interest  with  a  wide  variety  of  applications  ranging  from  communications1  to  laser  driven 
particle  accelerators.2'3  For  example,  recent  plasma  simulation  studies4’5  suggest  two  pos¬ 
sible  methods  by  which  the  frequency  of  an  electromagnetic  (EM)  wave  may  be  upshifted. 
In  the  first  method,4  the  plasma  density  through  which  the  EM  wave  is  propagating  is 
suddenly  increased  in  time,  while  the  second  method5  utilizes  the  interaction  of  a  plasma 
wave  wakefield  (having  a  phase  velocity  near  the  speed  of  light)2’3,6,7  with  a  short  EM 
pulse.8  In  this  Letter,  a  novel  method9  for  upshifting  the  frequency  of  a  laser  pulse  is 
proposed  and  analyzed,  which  utilizes  the  interaction  of  a  laser  pulse  with  a  relativistic 
ionization  front.  In  such  a  scheme,  an  intense  driving  electron  beam  (or,  alternatively,  an 
ionizing  laser  beam)  is  used  to  generate  an  ionization  front,  propagating  near  the  speed 
of  light,  with  the  plasma  density  rising  from  zero  to  some  high  value  n0  over  a  relatively 
short  distance  L.  A  “test”  laser  pulse  of  length  <  L  which  is  properly  phased,  such  that 
it  co-propagates  with  the  ionization  front  in  the  region  of  the  large  negative  local  density 
gradient,  will  be  continuously  frequency  upshifted  as  it  propagates.  This  mechanism  may 
provide  a  practical  and  efficient  method  for  tuning  the  frequency  of  a  laser  pulse,  in  which 
the  upshift  is  adjusted  by  varying  the  ionization  density  and/or  the  interaction  distance. 

The  physical  mechanism  by  which  an  ionization  front  may  be  used  to  upshift  the 
frequency  of  a  laser  pulse  may  be  understood  by  the  following.  Consider  a  laser  pulse  of 
frequency  w  and  wavenumber  k,  co-propagating  in  the  region  of  the  ionization  front,  both 
of  which  are  assumed  to  be  moving  at  approximately  the  speed  of  light.  Furthermore, 
assume  that  in  the  region  of  the  laser  pulse,  the  ionization  front  may  be  approximated  by 
a  linear  density  gradient  of  the  form 

r  ft,  £  <  -L\ 

n  ~  <  -nf,/L,  -L  <(,  <  0;  (1) 

U,  i  >  0; 

where  (  =  z  —  ct  and  L  is  the  length  of  the  ionization  front  (see  Fig.  1).  For  the  present 
discussion,  ft  =  no  and  L  =  Lo  are  assumed  to  be  constant.  Consider  a  laser  pulse 
centered  at  £c  =  — L0/2,  with  a  pulse  length  of  £,  where  L0  >  £  >>  X  =  2n/k,  and  of 
intial  frequency  u>n  >>  u>„n  =  (4ne2n0/mey^2 .  If  the  frequency  shift  induced  in  the  pulse 
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remains  small  compared  to  the  local  phase  velocity  of  the  pulse  (via  the  dispersion 
relation  u>2  ==  u2  +  c2k 2,  where  u>2  =  e2n/me)  is  approximately  vp/c  ~  1  +  u;p0n/2u>lno. 
Local  variations  in  the  plasma  density  n(£)  lead  to  variations  in  the  local  phase  velocity 
vp(£)  of  the  laser  pulse.  For  example,  the  local  phase  velocity  near  the  front  of  the  laser 
pulse  vp(£+)  will  be  less  than  that  near  the  back  of  the  pulse  vp(£_)  provided  n(£+)  < 
n(£_).  Hence,  the  individual  phase  peaks  in  the  pulse  ~  exp(tfcz  —  iu>f)  may  move  relative 
to  one  another  (i.e.,  closer  together  for  the  present  example).  The  resulting  change  in  the 
radiation  wavelength  may  be  estimated  by  A(r)  ~  A0  +  A vpr,  where  Avp  ~  A 0dvp/d(  is 
the  difference  in  phase  velocity  between  adjacent  phare  peaks.  A0  —  27rc/u>0  is  the  initial 
wavelength  and  cr  is  the  laser/ionization  front  propagation  distance.  This  gives  A/A0  — 
1— ct<j>20/2Lqu>1,  which  corresponds  to  an  increase  in  frequency  of  u>/u>o  ~  1+ctu>20/2L0w%. 
By  a  similar  mechanism,  a  plasma  wave  may  be  used  to  induce  shifts  in  the  frequency  of 
a  laser  pulse. s'8,9 

One  expects  the  above  estimates  to  hold  in  the  limit  in  which  the  frequency  shifts  re¬ 
main  small,  Au>  =  ui  —  wo  <<  o;0.  However,  it  is  possible  to  analytically  obtain  expressions 
for  large  frequency  shifts  (Au  >  u>)  in  the  limit  u>2/a;2  <<  1  via  the  wave  equation.  The 
transverse  plasma  response  current,  in  the  fluid  approximation,  is  J±  ~  —  envx,  where 
uj_  is  the  transverse  electron  fluid  velocity.  In  ID,  conservation  of  canonical  transverse 
momentum  implies  u_l  ~  ca,  which  is  the  quiver  motion  of  the  electrons  in  the  laser  field 
assuming  |a|2  <<  1.  Here,  a  =  eAx./mec2  is  the  normalized  vector  potential  of  the  laser. 
Using  the  variables  f,  ~  z  —  ci  and  r  —  t,  the  transverse  wave  equation  is  given  by 


\cd(  c2  dr )  dr  c2  no  h 


where  n(£,r)  is  the  plasma  electron  density  of  the  ionization  front  which  is  assumed  to  be 
known  and  independent  of  the  vector  potential  of  the  short  laser  pulse. 

The  wave  equation  may  be  solved  analytically  for  the  case  of  a  square  laser  pulse 
of  the  form  a  =  a  exp(iu;£/c),  for  —  t  <  t  <  0,  and  a  =  0,  otherwise.  Also,  a  linear 
ionization  density  profile  will  be  assumed,  of  the  form  given  by  Eq.  (1).  In  general,  the 
pulse  amplitude  a,  frequency  u;,  pulse  length  E,  ionization  density  ft  and  ionization  front 


length  L  are  all  functions  of  r.  In  the  limit  u;2/u;2  <<  1,  the  second  order  r  derivative  in 
Eq.  (2)  may  be  neglected,  and  one  finds  a  =  a0 u>0/u>  for  —  £  <  £  <  0,  where 


with  n0  =  n(r  =  0),  a0  =  a(r  -  0)  and  u>0  =  u>(r  =  0).  Notice  that  as  the  laser 
pulse  evolves,  the  laser  power  P  ~  \da/d(,\2  is  approximately  constant  in  t,  i.e.,  frequency 
upshifts  correspond  to  decreases  in  the  laser  amplitude,  a(r)  ~  l/w(r). 

By  integrating  the  wave  equation  in  £  over  the  longitudinal  extent  of  the  laser  pulse, 
an  equation  for  the  total  normalized  pulse  energy  U  =  J  d^\daj  d^\2  may  be  derived.  In  the 
limit  u>p/u>2  <<  1,  and  for  the  linear  density  gradient  given  by  Eq.  (1),  the  pulse  energy 
evolves  according  to  dU/dr  =  au;20n/2n0ci,  where  a  =  J  d£\a\2  >  0.  Hence,  the  laser 
pulse  gains  energy  from  the  ionization  front  as  it  propagates.  For  the  case  of  a  square 
pulse  of  length  £,  one  finds  £(r)/£0  =  cj(t)/lv0,  where  £0  =  £{t  =  0),  which  implies  that 
the  pulse  length  increases  as  the  frequency  is  upshifted. 

For  a  “fixed”,  externally  generated  ionization  front  (i.e.,  generated  by  an  ionizing 
pump  laser  or  electron  beam)  in  which  n  =  no  and  L  ~  Lq  are  independent  of  t,  Eq.  (3) 
gives 

w(t)  =  w0  (l  +  cru>20/L0Wo)1/Z  .  (4) 

Equation  (4)  indicates  that  asymptotically,  for  cr  >>  LoWq/u>20,  the  frequency  scales  as 
w  ~  upo{ct/ L0)1/2 ,  where  cr  represents  the  distance  the  pulse  has  propagated  through 
the  plasma.  Notice  that  frequency  shifts  Aw  ~  u>0  require  cr  ~  L0ul/u>20.  In  principle, 
Eq.  (4)  indicates  that  there  is  no  upper  limit  to  how  far  the  frequency  may  be  upshifted 
(assuming  that  the  ionization  front  and  laser  pulse  may  be  sustained  over  a  sufficiently 
large  distance  within  the  plasma),  however,  this  is  a  result  of  assuming  u0  ~  vg  ~  c,  where 
v0  is  the  velocity  of  the  ionization  front  and  vg  is  the  group  velocity  of  the  laser  pulse.  In 
practice,  u0  /  vg  which  implies  that  the  laser  pulse  will  “phase  slip”  in  £  out  of  the  region 
of  the  ionization  front  (-L0  <  (  <  0)  and,  thus,  will  no  longer  be  frequency  upshifted. 
Consider  a  laser  pulse  with  a  pulse  center  initially  located  at  £C(T  —  0)  =  £c0i  where 
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—  Lo  <  £cO  <  0-  Phase  slippage  will  cause  £c  to  evolve  in  r  according  to  d£c/dr  =  vg  -  v0. 
Assuming  v0  —  c,  vgjc  ~  1  —  iv20n(£c)/2u}2n0  and  n(£c)  =  -n0£c/L0,  where  u  =  o>(r)  is 
given  by  Eq.  (4),  gives  (c(t)/(cO  =  u{t)/co 0.  The  detuning  distance  crd,  defined  to  be  the 
propagation  distance  required  for  the  laser  pulse  to  phase  slip  out  of  the  ionization  front, 
(c(Td)  =  —  L0,  is  given  by  cTd  =  L0(e~2  —  ljcup/u;2,),  where  e  =  \£co/ L0\.  Inserting  this  into 
Eq.  (4)  gives  a  maximum  frequency  upshift  of  a>(rd)  =  u>0/e.  For  example,  a  pulse  initially 
centered  at  £ci)  =  L^jA  (e  =  1/4)  gives  a  final  frequency  of  u(Td)  =  4u>o  after  a  detuning 
distance  of  crd  =  15Lov£/u20. 

The  above  analysis  has  assumed  that  i)  the  velocity  of  the  ionization  front  is  constant, 
Vo  —  c,  as  may  be  the  case  for  an  ionization  front  produced  by  a  relativistic  electron  beam, 
and  ii)  the  density  profile  of  the  ionization  front  n(£)  was  assumed  independent  of  the 
propagation  time  r.  The  limit  placed  on  the  frequency  shift  as  a  result  of  phase  slippage 
may  be  circumvented  by  relaxing  either  of  these  two  conditions. 

If  the  ionization  front  is  directly  produced  by  the  laser  pulse  itself,  then  the  condition 
v0(r)  =  vg(r)  may  automatically  be  satisfied.  When  vq  =  vg,  there  is  no  slippage  and, 
hence,  the  frequency  of  an  intense  laser  pulse  may  be  “self-upshifted”  as  a  result  of  the 
interaction  of  the  pulse  with  its  self-generated  ionization  front.  For  a  square  pulse  which 
generates  an  ionization  front  with  a  uniform  gradient  dn/d£  =  -n0/Lo,  the  frequency 
evolves  according  to  Eq.  (4).  Since  «o  =  vg,  the  absence  of  slippage  implies  that  the 
upshifting  process  will  be  limited  by  some  other  mechanism,  such  as  laser  pulse  diffraction. 
As  an  example,  consider  a  100  fs  (£o  —  30  pm),  KrF  laser  pulse  (A0  =  2nc/w0  ~  0.25  pm) 
ionizing  hydrogen  gas  at  1  atm  (Apo  =  2tvc/u)po  —  4.5  pm).  In  order  to  double  the  laser 
frequency,  an  interaction  distance  of  cr  ~  3.0  cm  is  required.  Assuming  that  the  laser 
pulse  propagation  distance  is  limited  by  diffraction  implies  cr  ~  2Zr,  where  Zr  =  nr2t  / X0 
is  the  vacuum  Rayleigh  length  and  r,  is  the  laser  spot  size.  Setting  2 Zr  ~  3.0  cm  gives 
r,  ~  34  pm.  The  energy  required  to  ionize  this  volume  of  hydrogen  at  1  atm  is  about  1 1 
mJ.  Hence,  a  laser  pulse  energy  on  the  order  of  100  mj  should  be  sufficient.10  This  process 
of  frequency  self-upshifting  may  be  of  practical  significance,  however,  a  detailed  analysis  of 
this  process  requires  the  dynamics  of  laser  induced  ionization  be  included  seif-consistenti  v. 
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Alternatively,  the  limitations  resulting  from  phase  slippage,  which  occur  for  the  case 
of  an  externally  generated  ionization  front,  may  be  eliminated  by  gradually  increasing  the 
plasma  density  of  the  ionization  front  as  a  function  of  the  propagation  time  r.  By  correctly 
tailoring  the  ionization  density  as  a  function  of  r,  it  may  be  possible  to  control  the  laser 
group  velocity  such  that  v0  =  vg,  thus  eliminating  phase  slippage.  The  requirement  that 
vg  =  v0  is  constant  independent  of  r,  implies  that  density  must  be  tailored  such  that 
n(r)  ~  iu2(t),  i.e.,  Lotl{t)/ L{r)n o  =  u>2 /uig.  Solving  for  the  ratio  R  =  n(r)/T(r)  gives 
R/ Ro  =  exp(cru;p0/Lou>Q ),  where  i?0  =  R(t  =  0).  Hence,  the  ionization  profile  must  be 
tailored  such  that  hj L  is  an  exponentially  increasing  function  of  r,  which  may  be  achieved 
by  using  additional  lasers  to  control  the  ionization  rate.  For  such  a  density  profile,  Eq.  (3) 
indicates  that  the  frequency  is  exponentially  upshifted 

u’(r)  =  w0  exp(cru)20/2I0Wo).  (5) 

Zero  slippage  further  requires  that  the  ionization  front  velocity  be  held  constant  at  vq/c  — 
1  —  l&rolWpoy^Z/otUp,  which  may  be  satisfied  by  adjusting  the  energy  of  the  driving  electron 
beam. 

To  validate  the  analytical  results,  the  full  wave  equation,  Eq.  (2),  is  solved  numerically 
for  a  laser  pulse  having  a  half-sine  initial  profile,  |a|  =  |  sin(7r£/£0)|  for  -£0  <  (  <  0,  as 
shown  in  Fig.  1(a),  for  the  parameters  A0  =  2.5  pm,  Ap0  =  10  pm,  £0  =  25  pm,  L0  =  50  pm 
and  v0  =  c.  Equation  (2)  is  solved  for  two  cases:  (i)  h  =  n0  and  L  —  L0  are  constant,  which 
is  assumed  in  deriving  Eq.  (4),  and  (ii)  L  —  L0  is  constant  and  n  =  n0  exp(cTw^0 /2Lo<^o ), 
which  is  assumed  in  deriving  Eq.  (5).  The  results,  after  propagating  a  distance  cr  — 
(2  In  2)Z/0^po/^o  —  0-11  crn>  are  shown  in  Fig.  1(b)  for  case  (i)  and  in  Fig.  1(c)  for  case 
(ii).  The  results  show  excellent  agreement  with  the  analytic  predictions  (which  assumed  a 
square  pulse)  to  within  5%.  Specifically,  Fig.  1(b)  indicates  that  u>  ~  1.5u>0  (A  ~  0.65A0), 
|a|  ~  0.65|a0|  and  t  ~  1.5£o,  in  agreement  with  the  analytic  theory,  a  ~  1/w  and  t  ~  u>, 
where  u>  is  given  by  Eq.  (4).  Likewise,  Fig.  1(c)  indicates  that  u>  ~  2u>o  (A  ~  0.5A0), 
|a|  ~  0.5|ao|  and  l  ~  2£<j,  in  agreement  with  the  analytic  theory  where  uj  is  given  by  Eq. 
(5).  The  profile  distortion  and  the  “tail”  observed  behind  the  pulse  in  Figs.  1(b)  and  1(c) 
may  be  attributed  to  longitudinal  dispersion  arising  from  the  initial  half-sine  profile. 
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The  above  analysis  indicates  that  the  use  of  a  relativistic  ionization  front  to  upshift 
the  frequency  of  a  laser  pulse  may  be  of  practical  significance.  Physically,  frequency  up¬ 
shifts  arise  from  the  pulse  propagating  in  the  region  of  the  ionization  front  in  which  the 
local  density  gradient  is  negative.  This  produces  a  negative  gradient  in  the  local  phase 
velocity  within  the  pulse,  which  allows  the  laser  wavelength  to  decrease  and  the  frequency 
to  increase  as  the  pulse  propagates.  Interaction  with  an  ionization  front  having  a  uniform 
gradient,  dn/d£  =  —no/Lo,  evolves  the  frequency  as  described  by  Eq.  (4),  u>  ~  r1/2.  The 
laser  pulse  evolves  such  that  the  power  is  constant,  |a|  ~  1/cu,  whereas  the  pulse  energy  in¬ 
creases  due  to  increases  in  the  pulse  length,  i  ~  u>.  Furthermore,  phase  slippage  limitations 
may  be  overcome  by  appropriately  increasing  the  plasma  density  as  a  function  of  propaga¬ 
tion  time  t  which  may  lead  to  substantially  higher  frequencies,  i.e.,  as  described  by  Eq.  (5). 
Hence,  ionization  fronts,  by  varying  the  ionization  density  and/or  the  interaction  length, 
may  provide  a  convenient  method  for  continuously  upshifting  laser  pulse  frequencies.  In 
addition  to  this  mechanism,  an  ionization  front  may  also  be  used  as  a  relativistic  plasma 
mirror  upon  which  may  be  reflected  an  oppositely  directed  laser  pulse.9,11-13  The  reflected 
radiation  will  be  upshifted  by  the  relativistic  doppler  factor,  u,  ~  4-y where  u>,  is  the 
frequency  of  the  incident  radiation,  w,  is  the  frequency  of  the  scattered  (reflected)  radia¬ 
tion  and  =  1/(1  —  v\/c2)  is  the  relativistic  factor  associated  with  the  ionization  front. 
This  process  is  currently  being  investigated  and  will  be  the  subject  of  future  publication. 

This  work  was  supported  by  the  Office  of  Naval  Research,  the  Department  of  Energy 
and  the  Defense  Advanced  Research  Projects  Agency. 
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Fig.  I  —  The  initial  laser  pulse  and  ionization  density  are  shown  in  (a)  as  a  function  o f  (  =  Z  "  Ct.  Numerical  solutions 
of  laser  pulse  evolution  after  propagating  0.11  cm  are  shown  in  (b)  for  L  and  n  constant,  and  in  (c)  for  L  constant  and  n  an 
increasing  function  of  r. 
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